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SUMMARY 


This  paper  describes  shock-response  analyses  of  submerged  multiple  structures  with  two 
different  computational  methods.  Both  methods  approximate  the  presence  of  the  surround¬ 
ing  infinite  fluid  with  a  transmitting  boundary  based  on  the  first-order  Doubly  Asymptotic 
Approximation  [DAA\).  but  differ  in  their  treatment  of  multiple  scattering  by  the  struc¬ 
tures.  With  the  first  method,  the  DAA  boundary  is  placed  directly  on  the  structures'  wet 
surfaces,  which  is  strictly  valid  only  for  low-frequency  components  of  multiply-scattered 
waves.  The  second,  more  costly,  method  permits  a  more  valid  analysis  of  multiple  scatter¬ 
ing  through  finite-element  discretization  of  the  local  fluid  region.  Computational  results 
are  presented  for  simple  two-dimensional  problems  involving  two  circular  cylindrical  shells 
with  internal  masses.  The  results  produced  by  the  two  methods  are  often  in  close  agree¬ 
ment.  with  the  greatest  discrepancies  occurring  when  high-frequency  multiple  scattering 
is  important. 
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SECTION  1 
INTRODUCTION 


1.1  Background 

Analysis  of  the  response  of  underwater  structures  to  transient  excitation  has  long  posed 
a  formidable  challenge  to  engineers.  The  major  difficulty  in  solving  these  problems  is  the 
need  to  couple  structural  response  to  the  motion  of  the  surrounding  infinite  fluid,  which 
makes  conventional  finite-element  techniques  prohibitively  expensive. 

Early  efforts  concentrated  on  developing  analytical  solutions  for  a  single,  infinite  cylindri¬ 
cal  or  spherical  shell  see.  e  g..  Geers  ( 1975) ■ .  Although  exact  solutions  were  obtained, 
the  techniques  used  were  applicable  only  to  a  narrow  class  of  problems.  The  limitations 
of  analytical  approaches  prompted  the  search  for  general  purpose  numerical  tools.  An 
important  advance  in  this  direction  was  the  development  of  the  Doubly  Asymptotic  Ap¬ 
proximation  (DA  A)  for  treating  the  interaction  of  the  structure  with  the  surrounding  fluid 
Geers  (1971):  Geers  (1978):  Felippa  (1980)1.  In  essence,  the  DAA  replaces  the  infinite 
fluid  with  a  boundary  on  the  structure's  wet  surface  that  yields  correct  fluid  behavior  at 
high  and  low  frequencies,  and  effects  a  smooth  transition  between.  A  related  advance 
was  the  development  of  staggered  solution  procedures  Park,  et  ol.  (1977)  .  which  provide 
efficient  and  stable  means  for  integrating  separately  the  structure  and  DAA  equations  in 
time.  The  DAA  method  was  merged  with  a  staggered  solution  algorithm  to  produce  the 
USA  (Underwater  Shock  Analysis)  Code  DeRuntz.  et  al.  (1980)  .  which  has  been  found 
to  be  a  highly  efficient  and  versatile  too-. 

Unfortunately,  the  shock  analysis  of  multiple  structures  constitutes  an  application  where 
the  DAA  method  suffers  a  shortcoming,  for  the  DAA  does  not  properly  account  for  the 
high-frequency  components  of  multiply  scattered  waves,  although  the  low-frequency  com¬ 
ponents  are  treated  correctly.  This  limitation  is  overcome,  however,  by  application  of  the 
recently  developed  USA-STAGS-CFA  Code  Felippa  and  DeRuntz  (1984)  ,  which  places  a 
contained- fluid  field  between  the  structural  assembly  and  the  DAA  boundary.  The  con¬ 
tained  fluid  is  merely  a  local  portion  of  the  surrounding  fluid,  and  is  modelled  with  acoustic 
finite  elements.  Because  the  contained  fluid  extends  into  the  gaps  between  the  submerged 
structures,  it  is  possible  to  obtain  with  USA-STAGS-CFA  a  refined  analysis  of  wave  scat¬ 
tering  in  these  regions.  However,  a  high  price  must  be  paid  for  this  refinement.  In  addition 
to  the  burden  of  constructing  a  contained  fluid  mesh  and  integrating  the  contained-fluid 
equations,  the  time-integration  is  only  conditionally  stable,  which  limits  the  size  of  the 
integration  increment. 

1.2  This  Study 

In  view  of  the  increased  modelling  and  computational  demands  of  the  USA-STAGS-CFA 
code,  it  is  natural  to  inquire  if  adequate  solutions  for  multiple-structure  problems  can 
be  obtained  using  the  simpler  USA-STAGS  Code  DeRuntz  and  Brogan  (1980)  .  as  the 
theory  underlying  the  latter  is  valid  for  low-frequency  response.  The  investigation  of  this 
question  is  the  purpose  of  the  study  described  herein.  The  study  has  involved  comparing 


computational  results  obtained  with  both  the  USA-STAGS  and  USA-STAGS-CFA  Codes 
from  two-dimensional  shock  analyses  of  simple  multiple  structures.  The  complexities  of 
these  analyses  axe  such  that,  to  our  knowledge,  no  analytical  solutions  for  them  have 
been  obtained.  Consequently,  the  more  refined  USA-STAGS-CFA  results  are  used  as 
benchmarks  against  which  the  USA-STAGS  solutions  are  judged. 

1.3  Succeeding  Sections 

The  next  section  describes  the  governing  equations  for  the  finite-element  /boundary-ele¬ 
ment  methods  studied,  and  also  discusses  the  computational  techniques  used  to  solve  the 
equations.  Section  3  describes  the  two-dimensional  analyses  used  to  evaluate  USA-STAGS 
treatment  of  multiple-structure  problems,  presents  the  results  obtained,  and  states  some 
conclusions  drawn  from  the  evaluation. 


SECTION  2 

GOVERNING  EQUATIONS 


2.1  Introduction 

The  solution  of  underwater  shock  problems  for  multiple  structures  involves  as  many  as 
three  coupled  fields:  a  structural  field,  a  contained-fiuid  field  of  finite  extent,  and  a  DAA 
boundary  that  approximates  the  behavior  of  the  surrounding  infinite  fluid.  This  section 
describes  the  governing  field  equations  and  solution  methods  that  form  the  bases  of  the 
USA-STAGS  and  USA-STAGS-CFA  Codes.  The  intent  is  to  provide  a  summary  of  the 
relevant  theory:  for  a  more  extensive  discussion,  the  reader  is  referred  to  the  cited  refer¬ 
ences. 


2.2  Structural  Equations  of  Motion 

The  structural  field  is  treated  with  the  STAGS  Code  Almroth.  et  a 1.  (1980)  .  which 

employs  the  familiar  set  of  finite-element  equations  see.  e.g..  Zienkiewicz  (1977) 

Ms  xs  -  Ca  x,  -  K,x5  -  f  (2.1) 

where  x,  is  the  vector  for  the  structural  degrees-of-freedom  (DOF's)  and  f  is  a  load  vector 
arising  from  external  applied  forces  and  from  pressure  exerted  by  the  contained  fluid.  The 
matrix  is  a  diagonal  lumped  mass  matrix,  The  stiffness  matrix  K,  is  given  by: 

Ks  =  [  B'DBdV  (2.2) 

Jv. 

where  1*  is  the  structural  volume  and  D  is  a  constitutive  matrix  relating  stress  and  strain. 
The  matrix  B  relates  strain  at  an  interior  point  to  nodal  displacements  and  is  given  by 
B  -  LX'.  where  L  is  the  strain-displacement  operator  and  X  is  a  finite-element  shape- 
function  matrix  that  approximates  the  physical  displacement  r»($.t)  at  spatial  point  £ 
as 


x.(i.t)  *  X'($)x.(n  (2.3) 

The  damping  matrix  C.  in  ; 2 . 1 .  is  based  on  the  Rayleigh  damping  relation  Cs  =  qM,  - 
to  and  j  are  constants  .  which  ensures  the  existence  of  classical  normal  modes  in  the 
damped  structure. 

2.3  Equations  for  the  Contained  Fluid 

The  contained  fluid  is  an  acoustic  medium  of  finite  extent  that  lies  between  the  structure 
and  the  DAA  boundary.  The  treatment  of  this  medium  described  below  follows  that  of 
Newton  (1980)  and  Felipp a  and  DeRuntz  (19S4i. 

For  inviscid.  irrotational  motion,  the  fluid  displacement  if  can  be  derived  from  a  displace¬ 
ment  potential  c  as 
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where  p  is  the  fluid  density.  The  pressure  p  is  given  by 


P  =  v  (2.5) 

and  the  governing  equation  of  motion  is  the  wave  equation 

1 1  =  c2V\-  (2.6) 

The  constant  c  is  the  speed  of  sound  in  the  fluid  and  is  obtained  from  e2  =  Kip,  where  k 
is  the  bulk  modulus  of  the  fluid. 

Equation  (2.6)  may  be  discretized  by  application  of  the  finite-element  approximation 

t  =  N1^  (2.7) 

where  N  is  a  shape-function  vector  and  'I'  is  a  vector  of  nodal  values  of  w  .  An  equation 
for  'P  may  be  derived  with  the  Bubnov-Galerkin  method,  which  requires  that  the  weighted 
average  of  the  residual  error  in  (2.6)  should  vanish,  i.e.. 

I  N(tr"  -  c2V2tr)dV  =  0  (2.8) 

Jvf 

Application  of  the  divergence  theorem  to  (2.8).  followed  by  insertion  of  the  finite-element 
approximation  (2.7).  leads  to 


Q'l'  -  c2H$  =  c2b  (2.9) 

where  Q  and  H  are  symmetric  matrices  given  by 

Q  =  /  NN'dV  (2.10) 

Jv, 

H  =  /  (VN)(VNi‘dV  (2.11) 

J  v, 

and  the  vector  b  is  a  boundary  forcing  term  to  be  described  momentarily. 

The  matrices  Q  and  H  are  analogs  of  the  mass  and  stiffness  matrices,  respectively,  that 
appear  in  the  structural  field  equations.  As  in. common  structural  analysis  practice,  the 
response  calculations  are  expedited  through  replacement  of  the  consistent  Q  matrix  defined 
in  (2.10)  with  a  lumped  diagonal  matrix  denoted  Q.  Trie  lumping  is  accomplished  by 
placing  the  row  sums  of  Q  on  the  diagonal  of  Q. 

The  boundary  forcing  wctor  in  (2.9)  is  given  by 


i 

k 

\ 

» 


dB 


(2.12) 


where  B/  is  the  fluid  boundary  and  n  is  the  unit  normal  to  the  boundary  taken  positive 
outward.  The  integral  in  (2.12)  is  readily  simplified  to  a  more  convenient  matrix  expression. 
First,  observe  that,  from  (2.4). 


dv 

dn 


~b  - 
-pif  ■  n 


(2.13) 


where  is  the  fluid  displacement  on  the  boundary.  In  a  manner  analogous  to  the  dis¬ 
cretization  of  the  governing  equation,  the  boundary  displacements  may  be  interpolated 
from  nodal  values  as 


£)  *  K*b<  (2-14) 

where  Nfe  is  a  shape-function  vector  and  is  a  vector  of  boundary  displacements.  The 
use  of  (2.14)  and  (2.13)  in  (2.12)  leads  to 

b  -  pL  f  xbf 

where 

L /  =  -  /  NNl 

J  Br 

in  which  T b  is  a  diagonal  matrix  of  direction  cosines. 

It  is  interesting  to  note  that,  for  nodes  on  a  symmetry  plane,  b  —  0.  To  see  this,  merely 
observe  that  the  symmetry  boundary  condition  is  simply  £j  •  n  =  0.  Hence,  from  (2.13) 
and  (2.12).  b  =  0  on  that  boundary. 

2.4  Equations  for  the  DAA  Boundary 

If  a  structure  is  immersed  in  a  fluid  of  finite  extent,  the  contained-fluid  field  just  described, 
in  conjunction  with  appropriate  boundary  conditions  suffice'  to  define  fiuid  behavior  com¬ 
pletely.  On  the  other  hand,  if  the  fiuid  domain  is  infinite .  tru  presence  of  a  contained-fluid 
field  alone  is  adequate  only  if  the  field  extends  out  far  enougr.  to  avoid  interference  from 
waves  reflected  at  the  outer  boundary.  Thi-  st,-~aiied  '  island'  approach  iias  the  advantage 
of  being  "exact’  .  but  is  prohibitively  expensive  tor  a.,  but  'hurt-time  calculations.  An 
alternative  method  is  to  surround  an  abbreviated  contained-fluid  mesh  with  a  transmitting 
boundary .  which  allows  scattered-wave  energy  to  pass  mu  of  tne  contained-fluid  field.  The 
error  incurred  through  use  of  a  transmitting  houndar>  is  generally  small  if  the  boundary 
is  well  formulated. 

Asymptotic  Behavior.  The  transmitting  boundary  used  in  this  sTudy  is  based  on  the  Doubly 
Asymptotic  Approximation  (DAA)  Geers  1 1 97 1 ,  Geers  ( 1 97 8 ) :  Felipp a  (1980)  .  which  is 


TbdB 


A  b 


(2.15) 


(2.16) 
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exact  in  the  limits  of  both  low-frequency  and  high-frequency  motions.  It  is  appropriate 
to  begin  the  derivation  of  the  DAA  equations  with  descriptions  of  fluid  behavior  at  these 
extremes.  In  the  development  that  follow's,  it  is  supposed  that  the  DAA  boundary  is 
subdivided  into  a  mesh  of  boundary  elements,  whose  behavior  is  referred  to  nodes  located 
at  the  centroids  of  the  elements. 

Now  the  fluid  motion  may  always  be  represented  as  the  sum  of  an  incident  wave  and  a 
scattered  wave.  This  may  be  expressed  in  computational  vector  form  for  pressures  and 
normal  fluid-particle  displacements  at  the  DAA  control  points  as 


_  -  rJV  nSC 

Pd  -  Pd  ~  Pd 

V  ,  _  _  vsc 

Xd  —  Xd  -*■  Xd 


(2.17) 

(2.18) 


where  the  superscripts  IN  and  SC  denote  incident  and  scattered  waves,  respectively.  As 
the  incident  wave  is  specified  ab  initio,  it  is  only  the  effects  of  the  scattered  wave  that  must 
be  approximated. 

At  high  frequencies,  pressure  and  velocity  are  related  by  the  plane  wave  approximation 
(PWA)  Mindlin  and  Bleich  (1953)  .  written  as 


P !C  =  PcxSdr 


(2.19) 


The  physical  basis  for  this  relation  is  that,  in  the  high-frequency  limit,  acoustic  wavelengths 
are  much  shorter  than  the  characteristic  response  length  of  the  boundary  motion,  so  each 
element  of  the  boundary  can  be  thought  of  as  a  flat  plate  radiating  plane  waves  outw'ard. 

At  low  frequencies,  the  virtual  mass  approximation  (VMA)  applies  Chenock  (1970)  . 
w-ritten  as 


Adpf c  =  M,<x; 


(2.20) 


where  Ad  is  a  diagonal  matrix  containing  the  areas  of  the  DAA  boundary  elements.  The 
symmetric  matrix  Md  is  the  fluid  mass  matrix  for  computing  the  kinetic  energy  for  ir- 
rotational  flow  of  an  incompressible  fluid  that  is  excited  by  motions  normal  to  the  DAA 
boundary  DeRuntz  and  Geers  (1978)  .  The  physical  rationale  for  (2.20)  is  that,  in  the  low- 
frequency  limit,  acoustic  wavelengths  are  so  much  ionger  than  the  characteristic  response 
length  of  the  boundary  motion  that  the  fluid  behaves  as  an  incompressible  medium. 

DAA  Equation.  The  complete  DAA  equation  may  be  written: 


Mdpf° 


pc  AdPdL  ~  PcMdx'd( 


It  is  readily  shown  through  frequency-domain  analysis  that  (2.21)  embodies  the  proper 
asymptotic  behavior.  Taking  the  Laplace  transform  of  (2.21).  one  obtains 


sMdpfc 


pc  Adpd~ 


s1  pcMdxd( 


(2.22) 


IV.vCv. 
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where  the  bar  denotes  transformed  quantities  and  a  is  the  transform  parameter.  At  high 
frequencies,  or  large  s.  (2.22)  becomes 

p  sdc  =  spcxsdc  (2.23) 

which  leads  to  the  PWA  upon  back-transformation.  At  low  frequencies,  (2.22)  reduces  to 

Adp£c  =  s2Mdxfc  (2.24) 

which  becomes  the  VMA  upon  back-transformation. 

Applicability  of  DA  A  Boundary.  It  was  mentioned  in  Section  1  that  use  of  the  DAA  bound¬ 
ary  is  questionable  when  it  is  placed  directly  on  the  wet  surfaces  of  multiple  structures.  At 
low  frequencies,  the  VMA  can  be  validly  applied  to  multiple  structures  simply  by  calculat¬ 
ing  a  fluid  mass  matrix  that  couples  all  the  structures  through  the  fluid:  this  extension  has 
been  implemented  in  the  USA-STAGS  code.  Unfortunately,  the  high-frequency  PWA  is 
not  applicable  to  multiple  structures  because  the  PWA  assumes  that  scattered  waves  ema¬ 
nating  from  the  DAA  boundary  radiate  out  to  infinity.  This  assumption  is  clearly  violated 
for  multiple  structures  because  waves  scattered  by  one  structure  can  impinge  on  another 
structure,  so  the  waves  may  not  be  purely  outgoing.  One  can  argue  heuristically.  how¬ 
ever.  that  USA-STAGS  should  give  reasonable  results  if  the  high-frequency  components  of 
multiply  scattered  waves  do  not  significantly  affect  response. 

In  the  USA-STAGS-CFA  approach,  multiple  scattering  is  acommodated  within  the  con¬ 
tained  fluid.  This  is  clearly  a  more  rigorous  treatment  than  that  of  USA-STAGS.  although 
the  overall  analysis  is  still  approximate  because  the  interaction  with  the  infinite  fluid  is 
modelled  with  a  DAA  boundary.  In  this  instance,  however,  all  conditions  for  the  applica¬ 
bility  of  the  DAA  may  be  satisfied  by  fashioning  the  fluid  mesh  so  that  the  DAA  boundary 
is  everywhere  non-concave. 

2.5  Boundary  Conditions 

The  complete  description  of  the  behavior  of  the  component  fields  requires  the  specification 
of  boundary  conditions  that  couple  field  responses.  These  boundary  conditions  ensure  that 
force  and  displacement  compatibility  is  satisfied  at  interfaces  between  fields. 

USA-STAGS  Method.  In  tins  method,  the  DAA  boundary  is  placed  directly  on  the  wet 
surface  of  tiie  r-tructure  ana  only  one  interface  is  present.  Force  coinpatility  requires  that 
the  structural  force  vector  be  equivalent  to  the  pressure-force  exerted  by  the  contained 
fluid,  t.e.. 

=  -GdsAdpd  (2.25) 

where  f,  is  the  vector  of  nodal  forces  exerted  on  the  structure's  wet  surface  and  Gj,  is 
a  transformation  matrix  relating  structural  and  DAA  forces.  Application  of  the  contra- 
grauience  principle  Geers  ami  Ruzicka  !  19S4 .  then  yields  the  displacement  compatibility 
equation 


Xc  -  G^,XS  (2.26) 

USA-STAGS-CFA  Method.  In  this  method,  the  DAA  boundary  is  placed  on  the  exterior 
surface  of  the  contained  fluid  so  that  two  interfaces  are  present,  the  structure-fluid  interface 
and  the  fiuid-DAA  interface.  On  the  former,  force  compatilitv  requires  that  the  structural 
force  vector  be  equivalent  to  the  pressure-force  exerted  by  the  fluid,  j.e., 

f,  =  -G/.A/p,  (2.27) 

where  Gy4  is  a  transformation  matrix  relating  structural  and  fluid  surface  forces.  A/  gives 
the  tributary  area  on  the  structural  wet  surface  for  each  fluid  node,  and  py  is  a  vector 
of  fluid  surface  pressures.  Displacement  compatibility  on  this  interface  may  be  enforced 
by  entering  the  structural  displacements  at  the  appropriate  locations  in  the  vector  Xy 
appearing  in  (2.15).  An  alternative  procedure,  however,  may  be  used  that  eliminates 
the  need  to  calculate  and  store  the  matrix  Ly  defined  in  (2.16).  This  merely  involves 
application  of  the  contragradience  principle  to  obtain  Geers  and  Ruzicka  (1984) 

b4=pAyG'/4xs  (2.28) 

where  b4  is  the  contribution  of  the  structure  fluid  interface  to  the  boundary  forcing  vector. 
On  the  fluid-DAA  interface,  the  force  compatibility  relation  is 

pd~Gfdpf  (2.29) 

where  Gm  i-s  a  transformation  matrix  relating  DAA  control-point  forces  to  fluid  surface 
forces.  An  analysis  similar  to  that  used  in  the  derivation  of  (2.28)  leads  to  the  displacement 
compatibility  relation 

bd  =  pAdGydXrf  (2.30) 

where  bd  is  the  contribution  of  the  DAA  interface  to  the  boundary  forcing  vector  and  Ad 
is  a  diagonal  matrix  giving  the  tributary  area  on  the  DAA  boundary  of  each  fluid  node. 


2.6  Assembled  Response  Equations 

Based  on  the  preceding  development,  dynamic  fluid-structure  response  may  be  calculated 
by  step-by-step  numerical  integration  of  ordinary  differential  equations  in  time.  The  form 
of  these  ordinary  differential  equations  is  determined  by  the  computational  approach  se¬ 
lected.  In  the  following,  a  USA-STAGS  form  and  a  USA-STAGS-CFA  form  are  presented 
that,  in  the  opinion  of  the  authors,  provide  the  greatest  insight  into  the  solution  processes. 
However,  a  variety  of  condi'ion*-  have  dictated  that  neither  of  these  forms  be  implement- 
edu-ed  in  the  UiA-STA-  iS  and  USA-STAGS-CFA  Codes.  The  computational  approaches 
actually  employed  in  the  codes  are  described  ir.  DeRunrz  and  Brogan  (1980)  and  Felipp a 
and  DeRuntz  (1984)  . 


USA-STAGS  Form.  Here,  two  coupled  sets  of  ordinary  differential  equations  are  obtained 
as  follows.  The  first  set  is  assembled  by  introducing  (2.25)  and  the  first  of  (2.18)  into  (2.1) 
to  obtain 

M,x,  -  C,xs  -  Ksxs  =  -  GdaAd(PdN  ~  PdC)  (2-31) 

The  second  set  is  obtained  by  introducing  (2.26)  and  the  second  of  (2.18)  into  (2.21)  to 
produce 

MdpSdC  -  pcAdpsdC  =  Pc  MdfG^x,  -  x'ds)  (2.32) 

These  two  sets  are  solved  simultaneously  by  step-by-step  numerical  integration  in  time  for 
the  response  vectors  xs  and  pfc.  The  data  transfer  between  STAGS  and  USA  in  solving 
the  two  sets  is  illustrated  in  Figure  1. 

USA-STAGS-CFA  Form.  Here,  three  coupled  sets  of  ordinary  differential  equations  are 
obtained  as  follows.  The  first  set  is  assembled  by  introducing  (2.27)  and  (2.5)  into  (2.1) 
to  obtain 

]S4gX9  Csx*  K ax $  =  ~  G  j a A  j 4* s  (2.33) 

where  denotes  that  part  of  ^  pertaining  to  nodes  on  the  structure-fluid  interface.  The 
second  set  is  obtained  by  partitioning  as 

i>‘  =  {*'  <H‘0  (2.34) 

where  and  'i/d  denote  the  parts  of  ’F  that  pertain  to  interior  node,  and  to  nodes  on 
the  fluid-DAA  interface,  respectively.  The  introduction  of  (2.28).  (2.30)  and  the  second  of 
(2.18)  into  (2.9)  then  yields 

I  A/G/Sx,  | 

Q'F  -  c2H^  =  pc 2  O  (2.35) 

l  A dG'fd(x*c  -  x1/)  j 

Finally,  the  third  seT  is  obtained  by  introducing  the  first  of  (2.18).  (2.29).  and  (2.5)  into 
(2.21)  and  integrating  the  resulting  equation  twice  (with  quiescent  initial  conditions)  to 
produce 

pcMdx$c  =  Mrf(G fdifd  ~  PdS)  -  pcAdlGfdi/d  -  p‘dx)  (2.36) 

where  each  asterisk  over  denotes  an  integration  in  time.  The  equations  sets  (2.33). 
(2.35)  and  (2.36)  are  solved  simultaneously  by  step-by-step  numerical  integration  in  time 
for  the  response  vectors  x,.  i/  and  x^c.  The  data  transfer  among  STAGS.  CFA  AND 
USA  in  solving  the  sets  is  illustrated  in  Figure  2. 


2.7  Computational  Procedures 

The  USA-STAGS  and  USA-STAGS-CFA  codes  employ  staggered  solution  procedures.  Stag¬ 
gered  schemes  have  the  advantage  of  dealing  with  coefficient  matrices  that  pertain  only 
to  the  individual  component  fields.  These  matrices  tend  to  be  much  more  manageable 
than  the  motley  matrices  usually  generated  when  coupled-field  equations  are  merged.  In 
addition,  staggering  permits  the  optimum  assignment  of  a  time-integration  algorithm  to 
each  equation  set.  Yet  another  advantage  of  staggering  is  that  it  allows  the  individual  field 
processors  to  be  isolated  in  separate  software  modules. 

The  computational  procedures  used  by  USA-STAGS  and  USA-STAGS-CFA  are  discussed 
in  detail  by  DeRuntz  and  Brogan  (1980)  and  Feltppa  and  DeRuntz  (1984)  ,  to  which  the 
interested  reader  is  referred.  We  are  content  here  to  state  that  the  USA-STAGS  time- 
integration  procedure  is  unconditionally  stable  with  respect  to  time  increment  for  linear 
problems,  but  that  the  USA-STAGS-CFA  procedure  is  only  conditionally  stable.  The 
stability  limit  for  the  latter  is  roughly  given  by 


At  <  0(1  c) 


(2.37) 


where  I  is  the  smallest  distance  between  contained  fluid  nodes.  As  (2.37)  is  generally  more 
stringent  than  the  limit  imposed  by  integration  accuracy.  USA-STAGS  will  run  successfully 
with  fewer  time  steps  than  the  number  required  by  USA-STAGS-CFA  for  the  same  shock- 
response  problem. 
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SECTION  3 

TWO-DIMENSIONAL  ANALYSES 

3.1  Overview 

The  USA-STAGS  and  USA-STAGS-CFA  methods  are  here  compared  on  the  basis  of  two- 
dimensional  analyses  of  two  identical  infinite  cylindrical  shell  units  separated  by  a  distance 
of  one-half  radius  (Figure  3).  Each  shell  unit  consists  of  an  internal  oscillator  connected  to 
a  sandwich  shell  by  many  uniform  and  uniformly  spaced  springs.  Parameters  varied  in  the 
analyses  are  the  bending  stiffness  of  the  shell  and  the  fixed-base  natural  frequency  of  the 
internal  oscillator.  The  excitation  consists  of  an  incident  step-wave  oriented  either  side-on 
or  end-on  to  the  shell-unit  pair. 

The  multiple-shell-unit  analyses  are  supplemented  by  analyses  of  corresponding  single 
units.  In  addition  to  the  USA-STAGS  and  USA-STAGS-CFA  analyses,  the  single-shell 
studies  incorporate  an  analytical  form  of  the  DAA  approach  based  on  the  decomposition 
of  shell  response  into  Fourier  harmonics  Geers  (1974)..  The  single-shell  studies  serve  two 
purposes:  they  guage  the  extent  of  fluid  coupling  in  the  multiple-shell  problems,  and  they 
indicate  the  level  of  error  incurred  through  discretization. 

The  discussion  in  Section  2.4  suggests  that,  as  oscillator  fixed-base  natural  frequency  rises, 
agreement  between  USA-STAGS  and  USA-STAGS-CFA  results  should  deteriorate.  Also, 
it  is  reasonable  to  expect  greater  disagreement  between  the  results  for  an  end-on  wave 
than  for  a  side-on  wave  because  of  shadowing  effects  present  in  the  former  case. 

3.2  Description  of  Idealized  Structures 

The  sheli  units  are  most  conveniently  described  in  terms  of  the  dimensionless  parameters 
shown  in  Figure  4.  The  parameter  f!  is  the  fixed-base  natural  frequency  of  the  spring-mass 
system  inside  the  shell.  The  parameter  is  the  square  root  of  the  ratio  of  the  bending 
stiffness  of  the  sandwich  shell  to  the  stiffness  of  a  uniform  shell  constructed  by  removing 
the  core  and  fusing  the  inner  and  outer  surface  layers  of  combined  thickness  d.  This 
parameter  is  needed  in  order  to  model  the  ring  stiffeners  that  characterize  actual  pressure 
hulls  Geers  (1969).  Geers  (1974)  .  The  material  of  the  inner  and  outer  shell  layers  is 
characterized  by  its  density  (p«).  Poisson's  ratio  (i/).  and  plate  velocity  (c,).  the  last  given 
by  c*  —  Ea  Pali  -  v2).  where  Ea  is  Young’s  modulus.  The  values  given  in  Figure  4  pertain 
to  steel. 

3.3  Discrete-Element  Models 

Structural  Models.  Each  shell  is  modelled  with  a  single  ring  of  STAGS  410  shell  elements, 
with  72  elements  to  a  full  ring.  Symmetry  constraints  are  applied  as  required  to  the  shell 
edges,  including  axial  constraints  to  simulate  plane-strain  conditions.  Problem  symmetry 
is  exploited  by  using  half-ring  models  in  the  single-shell  and  end-on  wave  studies  and  by 
using  a  singie  full-ring  model  in  the  side-on  studies. 

The  identical  oscillator  springs  are  modelled  with  STAGS  200  beam  elements  that  extend 
from  every  node  of  the  sheli  model  to  an  oscillator-mass  node  at  the  centroid  of  the  model. 
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The  axial  stiffness  of  the  beams  is  adjusted  to  produce  the  desired  oscillator  frequency, 
while  the  bending  rigidity  is  zero.  The  axial  stiffness  of  the  beams  on  on  the  symmetry 
diameter  of  the  half-ring  models  is.  of  course,  half  the  nominal  value. 

Contained- Fluid  Models.  Different  finite-element  grids  for  the  contained-fluid  field  are 
required  for  the  single-shell,  end-on-wave,  and  side-on-wave  studies  (Figures  5-7).  The 
element  used  is  an  eight-node  isoparametric  brick  with  a  trilinear  interpolation  scheme:  at 
present,  this  is  the  only  fluid  element  available  in  CFA.  The  grid  geometry  for  the  single¬ 
shell  problems  was  generated  automatically.  The  grids  for  the  two-shell  problems  were 
first  sketched  by  hand  and  then  entered  into  the  database  with  a  digitizer. 

DAA  Models.  Each  DAA  mesh  matched  the  the  corresponding  structural  or  contained-fluid 
mesh  at  the  interface.  This  is  shown  in  Figures  5-7. 

3.4  Analysis  Procedures 

USA-STAGS  Analysis.  The  first  step  here  is  the  execution  of  the  STAGSC-1  Code,  which 
generates  the  structural  mass  and  stiffness  matrices,  as  well  as  geometry  data  and  other 
information  needed  by  the  USA  Code.  The  remaining  steps  involve  modules  of  the  USA 
Code  The  FLUMAS  module  accepts  the  geometry  of  the  DAA  boundary  and  generates 
the  fluid  mass  matrix  and  the  DAA-structure  transformation  matrix  Gda.  This  is 
followed  by  execution  of  the  AUGMAT  module,  which  generates  several  auxiliary  matrices 
used  in  solving  the  response  equations.  Finally,  the  TIMINT  module  performs  the  time 
integration. 

VSA-STAGS-CFA  Analysts.  The  first  step  here  is  identical  to  that  for  a  USA-STAGS 
analysis.  The  USA-CFA  versions  of  the  FLUMAS.  AUGMAT.  and  TIMINT  modules  are 
then  executed  in  sequence.  These  modules  are  similar  in  function  and  input  to  their  USA 
counterparts,  the  major  difference  being  that  geometry  data  describing  the  contained-fluid 
mesh,  rather  than  the  structural  mesh,  must  be  prepared  and  placed  in  a  file  to  be  passed 
to  the  FLUMAS  module  as  input. 

Time  Integration.  The  lengthiest  calculations  are  the  time-integration  runs.  The  costs  of 
these  are  inversely  proportional  to  the  time  increment,  which  should  therefore  be  set  as 
large  as  possible  within  the  constraints  imposed  by  accuracy  and  stability.  For  the  USA- 
STAGS  analyses,  where  only  accuracy  considerations  apply,  an  increment  of  .045a  c  was 
selected  on  the  basis  of  previous  experience. 

For  the  L S A-STAGr-CFA  calculations.  th<  nigne-r  accuracy  i>  achieved  by  setting  the 
increment  as  close  as  possible  the  Courant  stability  limit,  which  is  about  ,08a  c  for  the 
fluid  grids  used.  However,  coupled-system  stability  requirements  reduced  the  increment 
to  045a  c  for  the  smgie-shell  studies  and  025a  e  for  the  two-shell  studies.  Another  inte¬ 
gration  parameter  in  USA-STAGS-CFA  calculations  is  the  numerical  damping  parameter, 
which  is  here  taken  as  unity  Felippa  and  DeRunrz  (1984)  . 

The  responses  selected  for  evaluation  and  display  are  radial  velocities  at  the  front  and  back 
of  each  shell.  All  calculations  were  done  on  a  VAX  11  780  computer. 


3.5  Response  Results 

Unstiffened  Shell.  No  Oscillator.  The  first  shell  unit  to  be  considered  is  an  unstiffened 
shell  ("»  =  1)  with  the  internal  oscillator  absent.  The  incident  wave  is  a  step  wave,  with 
only  end-on  attack  considered  in  the  two-shell  analysis. 

The  single-shell  results  axe  shown  in  Figure  8.  In  addition  to  the  discrete-element  and  an¬ 
alytical  DAA  results,  the  figure  shows  exact  solutions  obtained  using  the  residual  potential 
method  Geers  (1969),  (1971).  (1972),  (1974)  .  The  exact  and  approximate  results  are  in 
close  agreement. 

It  is  appropriate  to  discuss  two  features  of  the  velocity  histories  in  Figure  8.  First,  the 
asymptotic  translational  velocity  of  the  shell  considerably  exceeds  the  fluid-particle  velocity 
of  the  incident  wave,  which  is  given  by  pcV  ■  P  =  1.  where  P  is  the  magnitude  of  the 
incident-wavepressure.  This  occurs  because,  with  the  internal  oscillator  removed,  the  shell 
is  very  positively  buoyant  ;see,  e.g.,  Geers  (1969)  .  Second,  the  small  “jumps’"  in  the 
histories  are  produced  by  extensional  waves  propagating  around  the  shell  that  gradually 
lose  energy  by  radiating  “creeping  waves’"  out  into  the  fluid  see.  e.g..  Geers  (1972)  . 

The  two-shell  results  for  end-on  attack,  along  with  their  single-shell  counterparts,  are 
shown  in  Figures  9  and  10.  With  regard  to  the  forward  shell  (Figure  9).  the  presence  of 
the  rear  shell  has  little  effect  on  response  at  the  front,  but  significantly  affects  response 
at  the  back  The  nature  of  the  latter  effect  is  as  follows.  When  the  incident  wave  and 
the  scattered  wave  generated  by  the  forward  shell  reach  the  front  region  of  the  rear  shell, 
that  area  moves  rapidly  inward,  generating  a  rarefactive  scattered  wave  that  propagates 
back  toward  the  forward  shell.  When  this  rarefactive  wave  reaches  the  back  of  the  forward 
shell,  it  pulls  that  region  radially  outward,  which  increases  velocity  response  markedly 
above  tha*  occurring  at  the  back  of  the  single  shell.  USA-STAGS  satisfactorily  captures 
this  effect  as  it  produces  results  in  close  agreement  with  those  of  USA-STAGS-CFA. 

Figure  10  shows  calculated  response  histories  at  the  front  and  back  of  the  rear  shell, 
along  with  corresponding  single-shell  results.  In  consonance  with  the  results  of  Figure 
9.  the  forward  shell  has  little  effect  on  the  response  at  the  back  of  the  rear  shell,  but  it 
significantly  affects  the  response  at  the  front.  Here.  USA-STAGS  does  not  do  so  well, 
producing  a  velocity  history  at  the  front  that  is  markedly  more  abrupt  than  its  USA- 
STAGS-CFA  counterpart  and  slightly  overestimating  peak  re-ponse 

Stiffened  Shell.  Low-Frequency  Oscillator.  This  shell  unit  consists  of  a  stiffened  shell  (">  = 
10 1  containing  a  low-frequency  oscillator  with  fixed-base  natural  frequency  Uaic  =  0.2. 
Side-on  as  well  as  end-on  step-waves  are  applied. 

The  single -shell  results,  shown  in  Figure  11.  show  tha:  here,  as  with  the  unstiffened  shell, 
there  is  close  agreement  between  the  USA-STAGS  and  USA-STAGS-CFA  results.  The 
slow  oscillation  in  Figure  11  reflects  the  interaction,  through  the  springs,  of  the  oscillator 
mass  and  the  combined  mass  of  the  shell  and  entrained  fluid.  The  frequency  of  this  motion 
is  somewhat  higher  than  the  oscillator  fixed-base  natural  frequency,  as  befits  reduced-mass 
oscillation.  Because  the  shell  unit  considered  here  is  neutrally  buoyant,  the  oscillation 
takes  place  about  the  fluid-particle  velocity  of  the  incident  wave.  It  is  interesting  to  note 
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that,  early  in  the  motion,  the  response  histories  in  Figure  11  agree  closely  with  their 
counterparts  in  Figure  8.  This  results  from  the  rather  high  flexibility  of  the  shell  walls  (be 
they  stiffened  or  unstiffened)  and  the  softness  of  the  oscillator  springs,  which  means  that 
early  response  is  dominated  by  inertial  and  membrane  effects,  which  are  identical  in  the 
stiffened  and  unstiffened  shells. 

Figures  12  and  13  show  response  histories  for  the  two-shell  configuration  excited  by  an 
end-on  step-wave.  In  general,  the  earlier  comments  on  Figures  9  and  10  apply  equally  well 
here.  Once  again,  the  highest  levelof  shell  interaction  through  the  fluid,  and  the  greatest 
discrepancies  between  the  USA-STAGS  and  USA-STAGS-CFA  results  occur  where  shell 
regions  are  in  close  proximity. 

Velocity  histories  for  the  side-on  wave  are  shown  in  Figure  14.  As  expected,  the  results 
are  generally  intermediate  between  those  for  the  single-shell  and  end-on-wave  cases  with 
regard  to  the  extent  of  shell  interaction  and  agreement  between  the  USA-STAGS  and 
USA-STAGS-CFA  results. 

A  somewhat  disturbing  feature  of  the  stiffened-shell  response  histories  is  the  presence 
of  small-scale,  high-frequency  oscillations  that  become  noticeable  at  ctia  ss  20.  These 
oscillations  suggest  slowly  growing  numerical  instability,  which  did  not  yield  to  treatment 
during  the  present  effort.  Hence  it  constitutes  a  high-priority  item  for  future  work. 

Stiffened  Shell.  High  Frequency  Oscillator.  This  shell  unit  consists  of  a  stiffened  shell  (->  = 
10)  containing  a  high-frequency  oscillator  with  fixed-base  natural  frequency  fia  c  =  1.0. 
Here  too.  side-on  as  well  as  end-on  step-waves  are  applied. 

The  single-shell  results  (Figure  15)  reveal  a  slight  deterioration  in  agreement,  relative  to 
that  exhibited  in  Figure  11.  between  USA-STAGS  and  Modal-DAA  results  on  the  one  hand, 
and  USA-STAGS-CFA  results  on  the  other.  This  is  consistent  with  an  earlier  study  Geers 
(1974)  i  comparing  exact  and  DAA  solutions,  where  it  was  found  that  DAA-based  methods 
tend  to  exaggerate  radiation  damping.  Here,  the  USA-STAGS-CFA  results  benefit  from 
the  two  layers  of  contained-fluid  elements  (Figure  5).  which  reduce  excessive  damping  by 
moving  the  DAA  boundary  away  from  the  shell's  surface. 

The  results  for  the  end-on  wave  (Figures  16  and  17)  reveal  similar  deterioration  in  agree¬ 
ment  between  the  USA-STAGS  and  USA-STAGS-CFA  results.  As  in  previous  compar¬ 
isons.  the  best  agreement  occurs  wi  er<  -hell  interaction  is  least  important,  viz.,  at  the 
front  of  the  forward  shell  and  at  the  bacs  of  the  rear  shell.  At  the  other  two  locations, 
discrepancies  are  pronounced,  especially  a:  the  front  of  the  rear  shell  (Figure  17).  where 
USA-STAGS  fails  to  account  for  early-time  Ihigh-frequencv)  shadowing  by  the  forward 
shell. 

Discrepancies  between  USA-STAGS  and  USA-STAGS-CFA  results  are  also  apparent  in  the 
side-on  wave-results  (Figure  18).  although,  in  comparison  with  the  end-on-wave  results, 
they  are  modest.  As  expected,  agreement  is  best  at  the  front  of  each  shell,  and  poorer  at 
the  back  of  each  shell,  where  some  multiple-scattering  effects  manifest  themselves. 

An  interesting  feature  of  the  “back"  velocity  histories  in  Figures  15-18  is  a  single-period 
oscillation  appearing  in  the  USA-STAGS  and  Modal-DAA  results  during  the  interval 


1  -  ct  Q  <  3  or  3  <  ct  >a  <  5.  The  fact  that  it  takes  20  Modal-DAA  circumferential  har¬ 
monics  to  capture  this  oscillation  in  Figure  15  is  evidence  of  the  short-structural-wavelength 
nature  of  the  phenomenon.  The  absence  of  this  oscillation  from  the  USA-STAGS-CFA  re¬ 
sults  reveals  an  important  drawback  in  the  use  of  finite-element  grids  to  propagate  transient 
waves.  The  USA-STAGS  method,  on  the  other  hand,  does  not  suffer  from  this  drawback, 
because  the  interaction  of  the  structure  with  the  infinite  fluid  is  handled  entirely  with 
boundary  elements. 

3.6  Conclusions 

The  USA-STAGS  and  USA-STAGS-CFA  results  for  transiently  excited,  multiple,  sub¬ 
merged  structures  are  likely  to  agree  satisfactorily  when  response  is  dominated  by  low- 
frequency  components.  Agreement  is  less  likely,  however,  when  intermediate-  and  high- 
frequency  components  are  significant.  In  these  situations,  the  use  of  USA-STAGS-CFA 
appears  to  be  mandatory,  although  the  presence  of  important  high-frequency  components 
requires  refined  meshing,  w’hich  in  turn  incurs  high  computational  costs.  In  most  cases.it 
should  be  possible  to  assess  the  applicability  of  the  simpler  USA-STAGS  code  to  a  three- 
dimensional  problem  by  comparing  results  obtained  from  the  application  of  both  USA- 
STAGS  and  USA-STAGS-CFA  to  related,  two-dimensional  evaluation  problems.  Such  a 
procedure  is  strongly  advised  in  production  analyses. 
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Figure  1.  Data  transfer  in  USA-STAGS  code. 
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Figure  2.  Data  transfer  in  USA-STAGS-CFA  code. 
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